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Abstract 



CsCuCi3 is a quantum triangular antiferromagnet, ferromagnetically 
stacked, with an incommensurate (IC) structure due to a Dzyaloshinskii- 
Moriya interaction. Because of the classical degeneracy caused by the frus- 
tration, fluctuations in CsCuCbj have extraordinarily large effects, such as the 
phase transition in longitudinal magnetic field (normal to the planes, parallel 
to the IC wavenumber q) and the plateau in q in transverse field (perpendicu- 
lar to q) . We argue that fluctuations are responsible also for the new IC phase 
discovered in transverse field near the Neel temperature Tn, by T. Werner et 
al. [Solid State Commun. 102, 609 (1997)]. We develop and analyse the 
corresponding minimal Landau theory; the effects of fluctuations on the frus- 
tration are included phenomenologically, by means of a biquadratic term. The 
Landau theory gives two IC phases, one familiar from previous studies; prop- 
erties of the new IC phase, which occupies a pocket of the temperature-field 
phase diagram near Tn, agree qualitatively with those of the new phase found 
exp er iment ally. 
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I. INTRODUCTION 



Compounds based on the triangular antiferromagnet (TAFM), particularly members 
of the ABX 3 family, have provided a wealth of interesting behaviour and indeed many 
surprises, as recently reviewedB. Because of the extraordinarily large effects of quantum and 
thermal fluctuation; CsCuCl, (with Neel temperaturj T N = 10.65 K aud zero-temperature 
saturation field H§ pa 30 TL ranks among the most interesting of these compounds; the 
surprises began 20 years agoS and continue still, in theory as well as in experiment. 

The magnetic properties of CsCuCl 3 are due to the Cu 2+ ions (S = 1/2); these form a 
stacked triangular lattice, to a good approximation. The interaction in the planes is antifer- 
romagnetic and therefore frustrated. Normal to the planes (in the chain or c direction), a fer- 
romagnetic interaction competes with a Dzyaloshinskii-Moriya (DM) interaction!!, giving an 
incommensurate (IC) structure with wavenumber q = qz in the c direction. In more detail, 
the classical, zero-temperature structure in zero magnetic field is the three-sublattice, ±120° 
TAFM structure; the spins lie in the planes, rotating by pa 5.1° per plane!. Application of a 
magnetic field yields a variety of interesting phenomena related to the classical degeneracy 
of the TAFM; recall that the classical ground state of the TAFM is continuously (and also 
discretely) degenerate, even in a magnetic field, and also that thermal fluctuations§Hi in 
classical models and quantum fluctuations^ break the continuous degeneracy in the same 
way, both selecting for example the colinear structure at H pa H§/3. 

In longitudinal magnetic field (normal to the planes, parallel to the IC wavenumber q), 
the low-temperature magnetization is discontinuousoll3'E3 at H pa OAHs, due to a novel, 
fluctuation-induced phase transitionBH: the umbrella structure is optimal at small H (due 
to a small, easy-plane anisotropy in the interplane exchang e§) and a coplanar structure 
is optimal at larger H (dueto quantum fluctuations). Other experimentJ^li support 
the Nikuni-Shiba analysist-ai^. In summary, CsCuCLj in longitudinal field appears to be 
understood, except that the transition at high temperature T has puzzling featuresB 

A transverse field (in the planes, perpendicular to q) gives more surprises. The be- 
haviour at low fields is conventional: q decreases quadratically with H, and the curvature 
increasesEl'0 with T. Classical (mean-field) theory@@ gives both of these and also a transi- 
tion to the commensurate (C) phase at if pa 0.47i?s ( a t T = 0); this transition has recently 
been observedil, but at larger H (pa 0.58i?s)- Classical theory cannot, however, explain any 
of the following: At low T, unusual behaviour occurs for fields near H = Hs/3, well below 
the IC— >C transition; structure is found in the mag netizatioiH, in the 133 Cs NMR shifti, 
in the IC wavenumber (Ref. ^ finds a plateau), and in ESR measurement sS. At high T, a 
second IC phase appears and increases with field0. 

At low T, the structure near H$/3 seems related to changes, induced by quantum fluc- 
tuations, in the structure of the TAFM near the same fieldu'EJ. Linear spin-wave (LSW) 
theory, which adds the leading quantum correction to classical theory, gives a plateau in 
the mae netization of the C state@ (as for the TAFM), a promising start. But there is 
another surprise, this time in theory. Not only does LSW theory of the IC phase fail to 
explain any of the other results, it actually provides a worse description of the IC phase 
than does classical theory, by predicting, a premature IC^C transitionil at H pa 0.32i?s- 
An innovativeDhenomenological theoryca of quantum fluctuations explains the existence of 
the plateauEHiJ in q, and also its level. The theoretical value {H pa 0.44if s ) for the field at 
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the IC^C transition!] is however too small, and the magnetization is not predicted well; 
likely the phenomenological theory can be improved. 

Both new fmdingsE-l near Tn, namely the increase of Tn with field (as in the TAFM@i@) 
and the second IC phase, are likely due to thermal fluctuations; neither has yet been treated 
theoretically. Because a microscopic or numerical treatment of fluctuations is out of the 
question for a vector-spin system with a nonsinusoidal IC structure, we use phenomenology. 
To treat a particular aspect of fluctuations near Tn, namely effects related to the breaking 
of the classical degeneracy of the TAFM, we add to the standard Landau theory a term 
biquadratic in the order parameter, as in our treatment^ of quantum fluctuations at T = 0; 
this term appears neither in the Hamiltonian nor in mean-field theoryl! at any T. Of course 
this term is not intended to include fluctuation effects in general (such phenomenology cannot 
possibly explain the reduction of Tn from the mean-field value of 35.5 K to the experimental 
value of 10.65 K). 

In qualitative agreement with experiment^, the Landau theory predicts a second IC 
phase to exist in a small T-H region near Tn. As discussed in the preceding articled and in 
Section IV, the theory explains other properties of the new phase, qualitatively; as discussed 
in the Appendix, however, it does not explain the increase of Tn with H. The new phase 
is the high-T version of a state which arose in the classical theory^; the state exists at all 
T in classical (mean-field) theory@@, but is never optimal. We argue that the new phase 
owes its existence to fluctuations; these are strong enough to overcome an energy difference 
in classical theory, just as quantum fluctuations in longitudinal field overcome the small 
anisotropy0£!. A unifying feature is that the two most striking ofthe experimental results 
in transverse field, namely the plateau@@ in q and the new phased near Tn, are explained 
using the same phenomenological treatment of fluctuations, in Ref. 29 and here respectively. 

Remarkably then, CsCuCl 3 displays a fluctuation-induced phase transition in transverse 
field, and a different fluctuation-induced phase transition in longitudinal fiekffi 



II. HAMILTONIAN 

The main interactions are described by the Hamiltonian 
'H = ^ [ — 2 J Sj n ■ Sj jn+ i — D z • (Sj n x Sj jn+ i) + J\ Si n ' Sfc n Q fJiB H x • Sj„ , (1) 



k 



where Sj n is the spin operator at the ith site in the nth a-b plane, z and x are unit vectors 
in the c and a directions, and the k sum is over the six, in-plane, nearest neighbours of the 
site in. The first term (oc Jo) is the isotropic, ferromagnetic exchange interaction between 
spins in nearest-neighbour planes, the second (oc D) is the interplane DM interaction, the 
third (oc J\) is the isotropic, frustrated, antiferromagnetic exchange interaction between 
nearest-neighbour spins in the a-b planes, and the fourth is the Zeeman energy in a field H 
transverse to the chains. We omit the easy-plane anisotropyS in the interplane interaction, 
the dipole-dipole interaction, and several other effects. The coefficients have been estimated 
previously&fiiffl; we use Jo = 28 K, J\ = 4.9 K and D = 5 K. The saturation field, above 
which each spin is aligned with the field at T = 0, is H s = 18JiS7(g/iB) ~ 30 T. 

At the classical level, the intrachain exchange term ( J ) favours states with spins parallel 
in adjacent a-b planes while the smaller DM term (D) favours states with spins in the planes 
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and rotating by 7r/2 per plane. At zero field, for all T < Tn, the spins lie in the planes, 
forming the 120° structure with three sublattices. The structure normal to the planes is 
helical; the wavenumber at H = is qo z where go — arctan (D/(2J )) « 27r/71. A transverse 
field deforms the helical structure, which becomes highly nonsinusoidal at higher H. 

For T = 0, the above Hamiltonian was investigated in the classical approximation!!, 
and the leading quantum correction was obtained using linear spin-wave (LSW) theory!. 
Neither theory can account for the structure observed near Hs/3, but a phenomenological 
treatmentl! of quantum fluctuations is largely successful. 

The extension of classical theory to T > (by mean-field theory) gives a phase diagram!! 
with only one IC phase, and so an understanding of the new IC phase near Tn seems 
to require including fluctuations at some level. A satisfactory microscopic treatment of 
fluctuations in CsCuCls is out of the question at T = 0, and they are even more difficult to 
treat for T > 0, leaving it seems only a jjhenomenological approach. For general T > 0, one 
could simply add the biquadratic termil to the mean-field expression^ for the free energy, 
but the coefficient would have to be adjusted as fluctuation effects increase with T, requiring 
a fit at each T or strong guidance from theory; this approach would be most reasonable near 
Tn, if the mean-field free energy were expanded to fourth order in the order parameter. We 
have chosen instead to use a fourth-order Landau theory. 



III. LANDAU THEORY 

The following describes a minimal Landau theory of CsCuCLj near Tn; usually, Lan- 
dau theory is reliable regarding the phase diagram, less reliable regarding the order of the 
transitions, and unreliable regarding fine details like the position dependence of the order 
parameter. We assume a structure with three sublattices in the a-b planes and with pe- 
riod L in the c direction (in units of the layer spacing); the restriction to integer L causes 
no difficultyii. We assume also that the spins remain in the a-b planes at all H and T. 
Curiously, the DM term is not sufficient for this!!!!; the easy-plane anisotropy0 helps of 
course, but fluctuation effects seem to be necessary!!. In mean-field theory, the free energy 
would be expressed in terms of the site-dependent magnetization {Sji), where j = 1,2,3 is 
the sublattice index and I — 1, . . . , L is the layer index. In Landau theory, the free energy 
is expanded in the order parameter m^, which is only proportional to (Sji). 

Explicitly, we use the following expression for the free energy F of the N spins, relative 
to the paramagnetic state at H = 0: 

yy 3 L 

F = TT J2 J2 [ \ «i m )i + «2 niji ■ m j+M - h x • + \ 71 - \ 7a (m j7 • m i+lj ,) 2 



31 j=i 1=1 



-Si mji ■ (m^j+i - mji) -8 2 z- (xaji x m j: i +1 ) . (2) 



This expression differs in major respects from the mean-field expansion to fourth order in 
m (or (S)). On the one hand, we omit from the expansion many terms of the same order 
as the ones we keep, second-order terms as well as fourth-order terms; on the other, and 
more importantly, we add a fourth-order term which does not appear in the expansion. 
All seven terms are essential. The third term in Equation (2) is the Zeeman energy; the 
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Landau parameter h is proportional to the magnetic field H. The terms with coefficients a\ 
and 71 are standard, while those with coefficients «2, 8\ and 82 result respectively from the 
in-plane antiferromagnetic interaction, the interplane ferromagnetic interaction and the DM 
interaction. The terms in ot\ and 8\ are adjusted so that the latter makes no contribution 
to the energy of the commensurate state. 

The remaining term, the biquadratic term — ^72 { m ji • m j+i,/) 2 ; appears neither in the 



Hamiltonian nor in the mean-field theoryc2l. It is introduced purely phenomenologically, 
to represent (to the extent possible) the effect of fluctuations (here thermal); the first such 
phenomenological use of the term was in Ref. In the TAFM, fluctuations (quantum or 
thermal) act to favour the colinear phase at intermediate fields, and therefore we require 
72 > 0; since the strength of the fluctuations is however unknown, 72 is otherwise adjustable 
(except that 72 < 0.5 for stability). In the TAFM, and in the C state of G3CUCI3, the 
fluctuations break the classical degeneracy, selecting one state from the many. Their effect 
in the IC phase of CsCuCls is far more subtle, for two reasons: first, the IC structure 
forbids the quantum selection possible in a C structure; second, the classical IC phase is 
unconventional, being well described as a continuous sequence of degenerate commensurate 
statesil. The fluctuations act to reorient the spins, giving a conventional IC phase with 
domain walls separating nearly commensurate regions. This is why LSW theory of the IC 
phase fails: a treatment beyond LSW theory is necessary to account for the reorientation. 

The above Landau theory is of course related to the mean-field theory of the ferro- 
magnetically stacked TAFM near Equation (3) of Ref. |36] (on the XY model), with 
Bi = 9Tn/5 for i = 1 to 6, reduces to our Equation (2) with the last three terms (coefficients 
72, 81 and 82) omitted, although some effort is needed to see this. As discussed in Ref. 



36, the colinear phase stable at intermediate fields cannot be obtained from this theory; it 



is obtained however if only B 4 is different from (less than) the other Bi, or if B 2 (or B 3 ) 
is different. Ref. |36] noted that some effects of fluctuations can be mimicked by allowing 
some Landau coefficients to depart from their mean-field values; these authors did not con- 
sider the biquadratic term explicitly, but their intention was like ours, to break the classical 
TAFM degeneracy. 

From the Hamiltonian, we expect the coefficients 02, 8\ and 82 to be proportional to 
6J1, 2J and D respectively. To reduce the number of parameters, we normalize the order 
parameter and the free energy so that «2 = 1 and 71 = 1. Then the interlayer coefficients 
are 81 = 2J /(6Ji) ~ 1.9 and 82 = Z?/(6Ji) ~ 0.17. Of the remaining parameters ai, h 
and 72, only a\ depends on T. At H — 0, the IC-paramagnetic transition at T N (0) occurs 
at oi\ ~ 1.015; cx\ = 1 is the upper limit of the Ciphase at H = 0. In Landau theory, 
the transition at 7n(0) is second-order. Experimented finds tricritical or weakly first-order 
behaviour; the latter is obtained in a recent Monte-Carlo analysis!^ of a related model. 

The unknown constants of proportionality involving m and h can be determined by 
comparing the Landau and mean-field theories. From Equation (2), the Landau energy (per 
spin) of the paramagnetic state at a% — 1 is |m 2 — hm + ^m 4 . The corresponding mean-field 
expression is found by setting f3 S 2 Jq = 1 in Equation (3.2) of Ref. 

Y c ( 18 J ^ s ) 2 - 2 ^b H(S) + ^^<s> 4 + • ■ ■) , 

where Jq = 4 J + 6 J\ and J = 2 J — 6 J\. Because the mean-field expression omits 
fluctuations, the term in 72 must be omitted from the Landau expression. On setting 
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(S) = am and comparing coefficients, one finds h = kH/Hs(0) where k = 3S/a and 
k 2 = 4 1 7 3 /(JqJ 1 ); the numerical value is k m 0.88. 

The Euler-Lagrange equations of the above Landau theory have many solutions, namely 
the paramagnetic (P) solution = mpx, several commensurate (C) solutions, and many 
incommensurate (IC) solutions; the following provides some background for the last. In 
classical theoryU, the magnitudes of the site magnetizations are fixed and the phases suffice 
for a complete description. A single IC solution is optimal at all fields below the IC— >C 
transition (predicted to occur at H « 0.47iJs); in this solution (called the 111 solution), the 
spins on all three sublattices wind through 2it over a period L of the IC structure. Many 
other IC solutions exist at fields below the transition. In the 110 solution, which is never 
optimal, the spins on only sublattices 1 and 2 wind through 2ir over one period, while the 
spins on sublattice 3 wobble about the field without winding. The 111 and 110 solutions, 
and many other solutions generated from them by forming composite solutions!!, become 
degenerate at the IC— >C transition, which is therefore a multi-phase point. In mean-field 
theoryil, for T > 0, the magnitudes of the magnetizations are no longer fixed; they can 
adjust to minimize the energy (for example by decreasing near a domain wall). The solution 
corresponding to the 111 solution is again optimal in all cases; it loses its winding character 
at larger T and H where the orbit in the m x -m y plane no longer encircles the origin. The 
infinite degeneracy at the IC^C transition remains at T > 0. 

The Landau-theory states corresponding to the 111 and 110 solutions are most easily 
described at H = (where the wavenumber is qo for both). For the first (ICi), 

vciji = mi [x cos(g ^ + fa) + y sin(g Z + fa)] for j = 1, 2, 3 (3) 

with fa = 0o + (j — l)27r/3. For the second (IC2), 

uiu = ~m 2 i = m 2 [x cos(g ^ + 0o) + y sin(g ^ + 0o)] , m 3l = ; (4) 

this is stable at H = only for ai > 0.99, but solutions at lower ai and H > are easily 
found. The first has the lower energy at H = 0, for all T. For the same amplitude, the 
first optimizes the a 2 term in the density of Equation (2), while the second optimizes the 
biquadratic term; the respective energies are — ^a 2 m 2 — |72m 4 and — a 2 m 2 — |7 2 m 4 . It is 
then possible that the second can have the lower energy, though orily for H > 0. 

For H > 0, some analytical results can be found at small field@'@, but full numerical so- 
lutions of the Euler-Lagrange equations are required at general values of H. Solutions were 
found by repeated linearization about trial solutions and solution of the linearized equations 
for the corrections. With increasing H, both solutions evolve, becoming increasingly nonsi- 
nusoidal; other solutions are found, but these are never optimal. In the second solution, the 
order parameter on the third sublattice increases from and wobbles about H with period 
L/2. Phase diagrams follow from comparison of the energies of the IC, P and C solutions. 

IV. PHASE DIAGRAM AND OTHER RESULTS 

Because Equation (2) cannot explain both the new IC phase and the increase of with 
H, and also because Equation (2) omits many terms of the same order as those kept, a 
detailed comparison with experiment is not attempted; fine details of the results should not 
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be taken seriously. The only parameter in the theory (apart from the unknown constant 
relating ot\ to T), is the coefficient 72 which determines the strength of the fluctuations. The 
value 72 = 0.2 used in the following was chosen, with guidance from experimentlll, to give a 
reasonable size to the IC2 region of the phase diagram, and also to the plateau in q/q for 
the ICi phase (as a function of h); 0.1 seems too small and 0.3 too large. Actually, the new 
phase appears even at 72 = 0, but only in a thin sliver (Ah < 0.024) of the phase diagram, 
with re-entrance and with no sizeable plateau. 

Figure 1 shows the theoretical phase diagram near Tn; again, a% is linear in T [cti = 1.015 
at T N (if = 0)] and h « 0.88H/H S (0). The two IC phases (ICi and IC 2 ), the P phase, and 
the C phase are optimal in the regions indicated; the C phase is the iri! = m 2 phase of the 
stacked TAFM. The new feature here is the IC 2 phase. 

The companion article!*] presents the strongest evidence for identifying the IC2 phase 
with the new IC phase discovered in Ref. |28|: the neutron-scattering intensity as a function 
of wavenumber is qualitatively that expected from the order parameter of the IC2 phase. The 
companion article compares theory and experiment in other respects as well. The following 
compares several aspects of the phase diagrams; it also presents results for the dependence 
of the IC wavenumbers on T and H, and for the order parameters in the two IC phases. 

Phase diagrams (qualitative aspects): In both Figure 1 and Figure 3 of Ref. [JT], the new 
phase appears in a small T-H region near Tn. Landau theory misses the increase of Tn with 
H and therefore also the nose; experimentally, the region of the new phase is shaped more 
like a croissant than a pocket. In both theory and experiment, the new phase appears only 
above some field, and a narrow tail extends to the low-T side. 

Phase diagrams (order of the transitions): The fields available in the high-T neutron 
measurementsll'0 were not sufficient to observe the C phase, nor of course the ICi — >C and 
C-P transitions; the IC^C transition was however observed at low T, at H ~ 18 T. 
In theory, all four of the other transitions in Figure 1 are first -order. ICi^IC 2 is strongly 
first-order, while the three transitions to the P state (ICi— >P near Tn, IC 2 — >P and ICi— >P 
at lower T) are weakly first-order (the free energies cross with almost the same slope); the 
ICi— >P transition is second-order at H — 0. 

In experiment, only ICj_-^IC2 is unambiguously first-order. The scanlll at 10.34 K can be 
interpreted in two waydS either the tail of the IC 2 phase was missed, or there is a first-order 
ICi— >P transition at H rs 12 T. The IC2— >P transition is almost certainly second-order, from 
the observation of critical scattering!!. 

Phase diagrams (other aspects): In theory, the IC 2 phase is not found at any T if h < 0.10 
or if h > 0.34, corresponding to H < 3.4 T and H > 11.6 T. In experiment, the lower limit 
is H = 4.3 ± 0.3 T; if the IC 2 phase does not appear at 10.34 K, then the upper limit is 
11.5 ± 0.5 T. The agreement is reasonable. Independent of the constant k relating h and 
H, the relative widths in the field variable are comparable: 11.5/4.3 ~ 0.34/0.10. The nose 
in the experimental phase diagram prevents a similar analysis for the temperature variable; 
for example, we cannot estimate reliably the upper T limit of the C phase (but none of our 
estimates disagrees with the data). 

Wavenumbers of the IC phases: Figure 2 shows theoretical results for the reduced 
wavenumber q/qo as a function of h, at four values of ot\. 

For «i = 1.005, there are two transitions as h increases, ICi— TC2 at h « 0.137 and IC2— >P 
at h « 0.202, both first-order in theory; at the first, q increases discontinuously to a value 
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less than the zero-field value go- I n theory, q/qo for the IC2 phase is roughly independent 
of h, for fixed ai- The dependence on a\ is stronger, but still weak; the value decreases 
roughly linearly with a 1; from q/q rs 0.98 at a\ = 1.01 to q/q ~ 0.88 at a\ = 0.995. The 
experimental value0 for q/q® in the IC2 phase is ~ 0.87; this is larger than in the ICi phase, 
as in theory; on the other hand, no dependence on T or H was observed. 
At lower T (ax = 0.98 and 0.96), q decreases as h increases, flattens out, bends over, and then 
drops discontinuously in a first-order transition to the P phase. At a\ = 0.93 (apparently 



corresponding to lower T than used in Ref. 31), q forms a reasonable plateau before rounding 



and falling to zero in a weakly first-order transition to the C phase; at slightly larger h, a 
second-order transition occurs to the P phase. 

Theoretically, the plateau in q/q occurs at ~ 0.6, almost independent of ax; the level of the 
plateau is reasonably robust (for 72 = 0.3, the plateau occurs at q/qo ~ 0.56). Figure 1 of 
the preceding article^ compares these results with the available data. Experiment finds a 
plateau (as in theory), at about the theoretical level (q/qo ~ 0.6), for both T = 10.34 K and 
T = 9.95 K; the latter data are slightly rounded at higher H, as in theory for cm < 0.98. 
Data were not obtained at low fields where q descends from the zero-field value to the 
plateau, preventing more detailed comparison with theory. 

Order parameters: Figure 3 shows the order parameter for the ICi phase, for Landau 
parameters a\ = 0.96 and h = 0.4 (near the middle of the plateau in Figure 2). At maximum 
m x , the configuration is almost colinear. Figure 4 shows the order parameter for the IC2 
phase, for a± = 1.005 and h = 0.18 (near the middle of the pocket in Figure 1). 
The preceding articled finds qualitative agreement between theory and experiment for both 
IC phases, but it is not possible to determine the order parameters uniquely from experiment. 
The theoretical order parameter in the ICi phase is too distorted (likely because of the loop 
at the right-hand side in Figure 3); that in the IC 2 phase agrees reasonably well, although 
the third sublattice is not visible in the available data. 



V. SUMMARY 

Landau theory with the biquadratic term (representing some effects of fluctuations) 
explains the appearance of the new IC phase foundHS near Tn. The new phase is the 
Landau-theory counterpart of the 110 state studied in Ref. |25], but stabilized by fluctu- 
ations. Only coarse adjustment of the only available parameter (the coefficient 72 of the 
biquadratic term) is needed to obtain qualitative (in some cases quantitative) agreement 
with experiment. 

In more detail: Landau theory finds a new IC phase to exist near Tn. It explains^ 
qualitatively the neutron-diffraction results in both IC phases. It does not explain the 
increase of Tn with field, but it explains other features of the phase diagram. It predicts 
moderately well the order of the transitions. It predicts that the wavenumber q of the IC2 
phase is larger than in the ICi phase, as observed; the experimental T dependence of q is 
however weaker than predicted. It predicts that the plateau in the wavenumber of the ICi 
phase occurs at q/q w 0.6, as observed. Theory and experiment agree qualitatively with 
respect to the order parameters. 
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APPENDIX A: EXTENDED LANDAU THEORY 

The Neel temperature Tn increases initially with field for the TAFM because of thermal 
fluctuations!!!. It increases also for CsCuCls (Refs. p8l , |31~D , for (one believes) the 



same 



reason. Since Equation (2) cannot explain the increase, it is natural to ask whether an 
extended Landau theory can do so. To investigate this question, we add the following 
fourth-order invariants^ to the square brackets in Equation (2): 



3 3 



k=l k=l 



+ \ 75 irijt ■ m j+1 j m fci ■ m k+h t + 7 6 (m,-^ ■ m j+1) i) ; (Al) 
fc=i 

our coefficients are related to those of Ref. |36| by 

Bi = 7i - 7;72 + 373 - 674 + §75 - 27 6 , 

B 2 = 71 - 272 + 47 6 , 
B 3 = 71 - 272 + 373 + I274 + 675 - 47e , 
B A = 71 - ±72 - 276 , 
B 5 = 7i + 72 + 373 + 374 - 375 + 7e , 

B e = 71 + 72 + 76 • (A2) 

For the stacked TAFM, at T sufficiently below Tn, the phase sequence with decreasing 
H must be: P phase — > C phase with mi = HI2 — > colinear C phase — > low-field C phase; 
a little analysis gives the requirements 72 > 0, -B3 > and B 5 > 0. If the P^C transition 
is second-order (as at h = 0, a\ — a 2 ), we find that the phase boundary is given by 



h 2 



oil ~ a l 

K 



B 

[ax + 2a 2 ) + 75-(«2 - cm 
B5 



(A3) 



Since B 5 > 0, the extended Landau theory cannot explain this effect of fluctuations (the 
increase of Tn with H) for the stacked TAFM; therefore we believe that it cannot explain 
either the increase of Tn with H for CsCuCl3. 

We also used the extended theory to determine several phase diagrams like Figure 1, for 
several different sets of parameters. The C, P and ICi phases always appear, as does the 
IC2 phase (unless of course 72 is sufficiently negative). Generally, the more complicated free 
energies (those with more of the ji parameters 7^ 0) give more complicated phase diagram 
(with for example re-entrant phases), but no new phases are found. 
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FIGURES 



FIG. 1. Phase diagram in the a\-h plane for 72 = 0.2. The Landau parameter a\ is linear in 
the temperature T, and h ~ 0.88H/Hs(0). The paramagnetic phase (P), the commensurate phase 
(C), and the two incommensurate phases (ICi and IC2) are optimal in the regions indicated. 

FIG. 2. Reduced wavenumber q/qo as a function of the Landau field h for four values of the 
Landau parameter a\. 

FIG. 3. ICi phase: orbit in the m x -m y plane for the order parameter on one of the three 
equivalent sublattices. The Landau parameters are ol\ = 0.96 and h = 0.4 (near the middle of the 
plateau in Figure 2); the period L is 116. The order parameters on the other two sublattices are 
displaced by I = ±L/3 from the first. 

FIG. 4. IC2 phase: orbits in the m x —m y plane for the order parameters on two of the three 
sublattices. The Landau parameters are a\ = 1.005 and h = 0.18 (near the middle of the pocket in 
Figure 1); the period L is 74. The outer loop is the orbit for one of the two equivalent sublattices; 
the orbit for the second is displaced by / = L/2. The order parameter on the third sublattice (inner 
loop) wobbles about the field with period L/2. 
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